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We show that the low-density strongly interacting electron liquid, interacting via the long-range 
Coulomb interaction, could develop a dispersion instability at a critical density associated with the 
approximate flattening of the quasiparticle energy dispersion. At the critical density the quasipar- 
ticle effective mass diverges at the Fermi surface, but the signature of this Fermi surface instability 
manifests itself away from the Fermi momentum at higher densities. For densities below the critical 
density the system is unstable since the quasiparticle velocity becomes negative. We show that one 
physical mechanism underlying the dispersion instability is the emission of soft plasmons by the 
quasiparticles. The dispersion instability occurs both in two and three dimensional electron liquids. 
We discuss the implications of the dispersion instability for experiments at low electron densities. 
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I. INTRODUCTION 

Recent theoretical worki on quasiparticle properties of 
interacting electron systems has indicated the possibility 
of a divergence of the zero-temperature quasiparticle ef- 
fective mass at a critical density, both in two-dimensional 
(2D) and three-dimensional (3D) quantum Coulomb sys- 
tems. In this paper, we investigate the nature of this 
mass divergence by considering the quasiparticle energy 
dispersion over a finite wavevector range around kp , ob- 
taining the interesting result that the effective mass di- 
vergence is not just a property of the Fermi surface (i.e. 
it does not happen only at kp), but is a dispersion in- 
stability where the whole quasiparticle energy dispersion 
around the Fermi surface becomes essentially almost flat 
- in fact, the instability first occurs at wavevectors far 
away from kp at densities well above the critical den- 
sity with the mass divergence eventually moving to the 
Fermi surface at the critical density. For densities lower 
than the critical density for mass divergence, we find that 
interaction effects drive the quasiparticle energy disper- 
sion 'concave' instead of the usual convex parabolic en- 
ergy dispersion of noninteracting free electrons, i.e., the 
energy actually decreases with increasing wavevector im- 
plying a 'negative' effective mass - the electrons have a 
'negative velocity' around the Fermi surface. Thus, the 
effective mass divergence within the ring diagram approx- 
imation for the electron self-energy reported in Ref. Q is 
actually a more general dispersion instability (i.e. the 
tendency of the quasiparticle dispersion of developing a 
flat band around kp) oi the type first discussed by Khodel 
and Shaginyan^ ~ the 'fermionic condensate' formation in 
the terminology of Ref. 0, fifteen years ago. 

We consider the standard jellium model for an elec- 
tron system with the electron-electron interaction being 
the usual '1/r' long-range Coulomb interaction (in 2D or 
3D) and the noninteracting kinetic energy being the usual 
parabolic dispersion with a bare mass 'to'. Such a system 
is characterized^ by a dimensionless interaction parame- 
ter, the so-called Wigner-Seitz radius = {■Kn)~^/^/aB 



(2D); (47rn/3)"i/VaB (3D), where n is the 2D or 3D 
electron density and gb = ft^/ (toc^) is the Bohr radius - 
Ts is both the effective interparticle separation measured 
in the units of Bohr radius and the ratio of the average 
Coulomb potential energy to the noninteracting kinetic 
energy. It was shown in Ref. Q that the on-shell quasi- 
particle effective mass, to*, diverges (both in 2D and 3D) 
when the quasiparticle mass renormalization is calculated 
within the leading-order approximation in the dynami- 
cally screened Coulomb interaction, or equivalently, in 
the infinite series of ring diagrams for the reducible po- 
larizability function. The critical r^-value (r*) for this 
effective mass divergence was found to be r* w 16 (2D); 
48 (3D). It was argued in Ref. that, although the spe- 
cific value r* is surely model-dependent (and one cannot 
expect the ring diagram approximation to give an exact 
or perhaps even an accurate value for the critical Tg), 
the qualitative fact that there is a quasiparticle effective 
mass divergence in strongly interacting 2D and 3D quan- 
tum Coulomb plasmas is a generic feature independent of 
the details of the approximation. In Ref. Q it was specu- 
lated that this effective mass divergence is the continuum 
analog of the Mott transition, or equivalently, a precur- 
sor to the Wigner crystallization (or perhaps a charge 
density wave instability). 



The other possibility is that this effective mass di- 
vergencei is the "fermionic condensation" discussed in 
Ref. 0, and fur ther elaborated, developed and discussed 
in Refs. I4ll5ll6ll7l l8l9l The "fermionic condensate" idea, pi- 
oneered by Khodel and collaborators, involves a fiatten- 
ing of the quasiparticle dispersion around kp at a, critical 
value of the interaction parameter, leading to interest- 
ing possibilities for the Fermi distribution function and 
the Fermi liquid theory. The inherent energy degeneracy 
associated with such a 'band flattening', where quasi- 
particles at different wavevectors have the same quasi- 
particle energy, could lead to various instabilities, and 
the 'renormalized degenerate' interacting system could 
in principle reorganize itself into a new (and perhaps an 
'exotic') phase such as a superconducting condensate^ 
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or a charge density waveiS. In such a situation involving 
a quantum phase transition to a non-Fermi hquid phase, 
the effective mass divergence is a signature of the emer- 
gence of the new collective phase, and the energy degener- 
acy associated with the band flattening (or equivalently, 
the plateau formation in the energy-momentum disper- 
sion) within the Fermi liquid phase itself is not a particu- 
larly critical issue (except as a precursor to the eventual 
quantum phase transition). The 'negative velocity' (i.e. 
decreasing energy with increasing momentum) that we 
find in our theory for > r* (i.e. below the critical 
density) obviously indicates an instability of the Fermi 
liquid, and our goal in this paper is to better understand 
this phenomenon. Note that the electron-electron inter- 
action we use in our model of 2D or 3D electron liquids 
is the realistic long-range Coulomb interaction. 

Since the effective mass divergence theoretically dis- 
covered in Ref. Q involves a very specific and ex- 
tremely well-understood many-body approximation^*ii, 
namely the ring-diagram approximation (sometimes re- 
ferred to as the RPA self-energy calculation or equiv- 
alently as the GW approximation^^) or its simple gen- 
eralizationsi^ incorporating approximate vertex correc- 
tions through local field corrections, we focus in the cur- 
rent work on a detailed calculation (within the same 
ring diagram approximation) of the full quasiparticle en- 
ergy dispersion i?(k) to understand the nature of the 
mass divergence by investigating whether the general- 
ized wavevector-dependent quasiparticle effective mass, 
m(k) = [{hk)~^dE{'k)/dk]~^, diverges only at k = kp 
or the quasiparticle dispersion is affected more drasti- 
cally for a finite range of k around kp akin to the band- 
flattening phenomenon leading to quasiparticle energy 
degeneracy. Wc find the latter situation to be the case 
with the whole energy dispersion being drastically af- 
fected by interaction effects at large Vg. We also find 
that within our leading-order self-energy calculation (in 
the dynamically screened Coulomb interaction) the dis- 
persion instability and the associated effective mass di- 
vergence arise, at least partially, from the emission of col- 
lective modes ("plasmons") by the quasiparticles, which 
could happen at rather low wavevectors {k > kp) in the 
strongly interacting regime of large (where the mass 
divergence phenomenon occurs). Our identification of 
a specific physical mechanism contributing partially to 
the mass divergence and the dispersion instability phe- 
nomena, namely the spontaneous emission of collective 
plasmon excitations by the quasiparticles leading to a 
negative quasiparticle velocity (or equivalently, a strong 
nonlinear "dip" in the quasiparticle dispersion) , suggests 
that the dispersion instability and the effective mass di- 
vergence is a generic feature of strongly interacting quan- 
tum Coulomb electron systems and is not a artifact of the 
ring diagram approximation. 

The rest of this paper is organized as follows. In ^ 
we describe the theory and formalism; in IIIII we present 
our results for the quasiparticle energy dispersion; and 
we conclude in IIVI giving a detailed discussion of the im- 



plications of our results. 



II. THEORY AND FORMALISM 

In this section we present the formalism we are going 
to use in our work. Our goal is to calculate the quasi- 
particle energy dispersion as a function of the interaction 
parameter . The central quantity needed for this calcu- 
lation is the electron self-energy function. Without any 
loss of generality we assume the electron system to be 
spinless for our calculations since spin plays no role in 
the theory as Coulomb interaction is spin independent. 
We choose h ^ 1 throughout, which makes wavevector 
and momentum (as well as energy and frequency) equiv- 
alent. 



A. Random Phase Approximation 

We examine the jellium electron system with long- 
range Coulomb interaction between electrons at zero 
temperature. Within random phase approximation 
(RPA)5iiiii2ii2ij the part of the ground state energy in- 
troduced by Coulomb interaction can be denoted by the 
Feynman diagrams shown in Fig. ^ Following Landau's 
approach, the quasiparticle energy can be obtained by 
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(1) 



where Uk is the distribution function at momentum k 
and Eg is the ground state energy of the system. 
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FIG. 1: The Feynman diagram for the Coulomb interaction 
contribution to the ground state energy within RPA. The cir- 
cles are polarization bubbles, the wiggly lines are the bare 
Coulomb interaction, and the solid lines the noninteracting 
electron Green's function. 

Note that the distribution function enters into the ex- 
pression of ground state energy Eq through the nonin- 
teracting Green's hmction 



G(")(k,c^) = 



1 - «k 



17] 



(2) 



where = fc^/(2m) — Ep is the non-interacting elec- 
tron energy with Ep as the Fermi energy or the non- 
interacting chemical potential, and 77 is an infinitesimal 
positive number. It is easy to obtain the variational 
derivative of the Green's function as 



G(") (q, uj) = 2TTid{uj - a)'5(k - q). 



(3) 
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Graphically, taking the rik variational derivative of a 
quantity simply means cutting one solid line of the Feyn- 
man diagram and taking the external momentum and 
frequency to be on-shell (i.e. uj = ^k)- We show the 
corresponding Feynman diagram for the self-energy in 
Fig. [3 We emphasize that the RPA or the ring-diagram 
approximation (which is appropriate for electron liquids 
interacting with the long-range Coulomb interaction) as 
shown in Fig. ^ necessarily implies that the on-shell self- 
energy approximation is used for calculating the quasi- 
particle energy dispersion (Fig. \^ since all energy and 
momenta in Fig. ^ correspond to the nonintcracting sys- 
tem. 



+ 




FIG. 2: The Feynman diagram for the self-energy within 
RPA. The thick wiggly line denotes the dynamically screened 
Coulomb interaction. 

The second order derivative of the total ground energy 
is referred to as Landau's interaction function: 



/(k,k') 



drikdrik 



(4) 



The Feynman diagram for the interaction function within 
RPA is shown in Fig. El 



where q = |k — k'| is momentum transfer between the 
two interacting electrons, and k is the appropriate inverse 
screening length. The minus sign in Eq. © reflects the 
exchange character of the Coulomb interaction between 
two fermions. Because the amplitude of faiq), given by 
Eq. ©, is maximal at g = 0, i.e. at 6* = 0, where 0036* > 
0, it produces a positive contribution to the right-hand 
side of Eq. (O and tends to decrease m* . Because the 
first term fa dominates over the other two terms fb and 
fc in Fig. 131 at small Ts , thus m* < to at small rg (see 
Fig. 6 of Ref.lnl). 

However, the sum of the ft and fc terms in Fig. Ohas 
the same sign as fa in Eq. lO, but its maximal amplitude 
is achieved at k = — k', i.e. at = tt, where cos 9 < (see 
Fig. 10 of Ref. ITU) . Thus, fb and fc produce a negative 
contribution to the right-hand side of Eq. ||SJ) and tend 
to increase m*. Because the strength of the and fc 
terms ( "two-plasmon exchange" ) increases relative to fa 
( "one-plasmon exchange" ) with the increase of , the 
effective mass m* starts to increase and becomes greater 
than TO. This effect was already seen in Fig. 6 of Ref. ini 
However, the numerical calculations of Ref. O stopped 
at a moderate ~ 4. In Ref. essentially the same 
calculations were extended to greater r^, and divergence 
of TO* was found. The divergence of to* occurs at the 
point where the interaction term in the right-hand side 
of Eq. ||SJ) becomes equal to the bare term. 

An alternative, but in some sense similar scenario of 
mass divergence was proposed in Ref. based on a phe- 
nomenological assumption that /(g) has maximal ampli- 
tude at g ^ 2kp, where cos 9 < 0, due to proximity to a 
density-wave instability. 
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FIG. 3: The Feynman diagram for the Landau interaction 
function /(k, k') within RPA. 

Landau's interaction function |0J determines renormal- 
ization of the effective mass to* relative to the bare mass 



1 

TO* 



C / f{9)cos9do, 



(5) 



where 9 is the angle between k and k', do the element 
of solid angle along k' in 3D and d9 in 2D, and C 
kp/i^irf in 3D and C = 1/{2t:)^ in 2D. 

The contribution (a) in Fig. El to Landau's interaction 
function is nothing but the static screened Coulomb in- 
teraction 



faiq) 



27re^ 

q + K 



in 2D, fa{q) 



in 3D, 



(6) 



B. Dielectric function 

The key quantity in evaluating the self-energy is the 
dynamically screened Coulomb interaction u(q, w) = 
Vq/e{q,Lu), where Vq is the Coulomb interaction in mo- 
mentum space and e(q, w) is the dynamical dielectric 
functiori2iiiii24i^. The Coulomb interaction is given by 
Vq = 27reVg (2D) and ine'^/q'^ (3D) with g = |q| be- 
ing the appropriate 2D or 3D momentum. The dielectric 
function is given by the infinite geometric series of the 
nonintcracting ring diagrams where each ring is just the 
nonintcracting electron polarizability. The 2D and 3D 
dielectric functions e{q,Lu) = 1 — VqIl{q,Lu) are given in 
Stern^^ and Linhard^^- respectively, where H(g, uj) is the 
polarizability, which can be denoted as one ring or bubble 
as in Fig. [Hand Fig. El 

In actual calculations it is conventional to express 
all the expressions in terms of the dimcnsionless units 
rg. The relation between rg and Fermi vector is rg — 
l/iakpaB), where a^^^^ = 1/2 and a'-^"^ = (97r/2)-i/3^ 
For completeness we write down the expression of 2D and 
3D e(g, oj) in the units of 2kp = AEp = 2m ~ 1, i.e. we 
use X — qj i^kp) to denote momentum and u — uj/(4:Ep) 
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to denote energy. For 2D u > case we have: 



C. Quasiparticle self-energy 



Re e(x,u) = 1 + — ^< 1 
2x 



The one-loop self-energy (Fig. ^ at zero temperature 
can be written as 



2x^ 



-— rSign(a::^ + u)e[{u + x^ f - x^]J{u + x^ f - x^ 
2x^ 



S(k,c^) 



di^q f dv Vq 



{2tt)'^ J 27rie(q,z/) 
•G(°)(k + q,L^ + i.). (9) 



Im e{x,u) = ^^e[x^-iu- x^fWx^-{u-x^)^ 



e(x, iu) = 1 + 



2x 



— Jx^ - x'^ -u + v/(x4 - x2 - u)2 + 4M2a;4( 7) 



where Q{x) = 1 when a; > and otherwise. For 3D 
M > case we have 



Re e{x, u) — 1 



■kAx'^ 



q^ + q 



q^ ~ q 
q^ + q 



q^ -q 



2x3 




FIG. 4: Contour in the frequency v plane. The cross denotes 
the possible position of the pole of Green's function. 



Due to the difficulty with the principal value integra- 
tion and singularities of l/e(k, a;) along the real axis in 
Eq. lO, it is advantageous to follow the standard proce- 
dure and deform the frequency integration from the real 
axis to the imaginary axis. After choosing the contour of 
frequency integration as in Fig. ^ we consider the inte- 
gration 



(10) 



Im e(x, u) 



e(a;, iu) 







ttAx^ 
1 , 



(x - > M > 0) 
u)^] {x^ + X > u> — a;|) 
(Otherwise) 



4a;3 



2 4 
X — X 



u^] In 



{q^ + q? + 



{q^ — q)"^ + 



u u 
arctan(— ^ ) — arctan(— ^ ) 



"q^ + q' 



q^ - q' 



where it;(q, v) is any complex function that is analytic 
in the upper and lower half of the complex plane which 
satisfies the condition that ?«(q, i') — > as \v\ oo. On 
the one hand, pU|l equals the integration of the integrand 
along real and imaginary axes since the integration along 
the curved part of the contour F (Fig. ^ vanishes. On 
the other hand, (|10|l equals the residue due to the pole 
of Green's function 

- [1 - nF(^q+k)]e(tJ - ^q+k)w(q, Cq+k - UJ - if]) 
+ '^i^(Cq+k)©(Cq+k-^^)w(q,^q+k-C^ + W7), (H) 



(k) where npix) is the Fermi distribution hmction (At T = 
J npix) = 1 for X < and otherwise). Now by setting 



In both 2D and 3D we have the relation e{x, — u) 
e*{x,u) (complex conjugate) and e{x,—iu) = e{x,iu). 



u'(k, Lo) 



(12) 
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we have 



d'^q [ dv , -.^ , , N 
TTTTd i T^— Wqw(q,i/)Go(q + k,j/ + 

(ZTTj" J Im 



f d'^q 



(13) 



• ©(^^ - Cq+k)w(q, ^q+k - CO - i?7) 
+ "-F('?q+k)0(Cq+k - C^)w(q, ^q+k - + i77)| 

r d'^q r dv w{q,w) 
J {2ttY J 27r iz/ + L^-eq ■ ^ ^ 



?q+k 



Note that the first term in Eq. usually named the 

exchange part, is a singular term when temperature is 
zero. Thus we rewrite the first and the second term so 
that they are well defined. Real and imaginary parts of 
the self-energy E(k,w) = S'(k,w) -|-iE"(k, w) can then 
be written as 



S'(k, ^) = - j -^v,Q{2mu: + 4 - |q - kp) 



d'^q 



TJoRe 



' e(q,^q-k-c.) 
Q{2muj + fc| - |q - kp) - e(fc| - |q - kp) 



(2^ 



2^^ 



e(q, 



1 



W + CJ - $q+k 



(15) 



Jm 



1 



e(2ma; + k} 



e(q,Cq-k -uj) 
q-kn-e(4-|q-k|il6) 



To be consistent with our leading-order one-loop ap- 
proximation in the dynamically screened Coulomb in- 
teraction (Figs, n and O we should calculate the self- 
energy only within the on-shell approximationi*ii instead 
of solving the full Dyson's equation. After putting w = f k 
(i.e. the on-shell approximation), we express the above 
equations in terms of r^, while using 2kF as the unit of 
wave- vector, and AEp as the energy unit. For 2D real 



part of self-energy Ek = 5](k, w = ^k) we obtain 

'-y 



AEf 



dx / dO- 



Ji/2 Jo \/x^ — 2xy cos{9) + 

■ Re ^ ^ 

e I x^ — 2xycos{6) + y^, x^ — y"^] 



V2Tr 



dx / du 



e(x, ixu) 



(x^ -V?- 4y2) + ^(x2 ~ V? - 4y2)2 _|_ 4y2y2 



where y = k/(2fci?). Similarly for 3D we have 



dx / dO- 



J 1/2 Jo x^ ~ 2xy cos{0) + y'^ 

■ Re ^ r- 

e I ^ x^ — 2xycos{6) + y^, x^ ~ y^j 



arc 



C30 /"OO 

dx / du 



47r^y Jo JO 
, ■(2y-a;)2+u2 
In 



e(x, ixu) 



(18) 



_{2y + .x)2 + u2 
The 2D imaginary self-energy can be written down as, 



dx / d0- 



AEf tt J Jo ^x^ — 2xycos{0) + y^ 

1 



•Im 



(19) 



^y^x^~^-^an7cos(^y^r^, x^ — y^^ 
and for 3D imaginary quasiparticle self-energy we have 



dx / de- 



AEf tt ii/a^'io "''a:2-2xzjcos(6')+y2 

1 



•Im 



e (^y^ x^ — 2xy cos(6') -I- _ y2~j 



(20) 



D. Quasiparticle energy dispersion and damping 

The quasiparticle dispersion is given simply by adding 
the on-shell real self-energy to the nonintcracting electron 
energy 



-E'k = ^k + Ejj. 



(21) 



The quasiparticle damping rate, which is non-zero away 
from the Fermi surface fc 7^ fc_F, is given by the on-shell 
imaginary part of the self-energy 



rk = |imE(k,ek)| = e;:. 



(22) 
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The renormalized quasiparticle effective mass is given by 



m*(k) 



1 dEk 
k dk 



(23) 



It is clear tliat we get m* = m, the constant bare electron 
mass, if we use the noninteracting energy = fc^/(2rn) 
for E^, i.e. if we put the self-energy to be zero. Putting 
k = kp one gets the quasiparticle effective mass at the 
Fermi surface, as calculated in Ref. 1'. For an arbitrary 
k ^ kp, m(k) defines the dispersing quasiparticle mass, 
i.e. the generalized wavevector-dependent quasiparticle 
effective mass. 



III. RESULTS 
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In this section we present our calculated results show- 
ing the quasiparticle energy dispersion for both 2D and 
3D electron systems with long range Coulomb interac- 
tions. The dispersion instability at large will be obvi- 
ous in these results. 

In Fig.|31we plot the calculated quasiparticle energy as 
a function of momentum for different Ts values in a 2D 
system. It is clear that the quasiparticle effective mass 
(m* — kF{dEk/dk)~^\k=kF) diverges at some low den- 
sity at critical value r* ~ 18. As increases, not only 
the energy dispersion of the quasiparticles around Fermi 
surface becomes flatter (i.e. relatively independent of 
momentum), but the dispersion in the whole momentum 
space is fundamentally changed. The quasiparticle en- 
ergy dispersion in a 3D electron system is presented in 
Fig. El showing the dispersion "flattening" for a range of 
wavevector fc at a critical r* ~ 50. 
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FIG. 5: (Color online) 2D quasiparticle energy as a function 
of momentum at values ranging from to 20. Inset: zoom- 
in of the spectrum with k/kp < 2. The dashed line indicate 
the Ts value at which the effective mass at the Fermi surface 
shows divergence (r* ~ 18). 



FIG. 6: (Color online) 3D quasiparticle energy as a function 
of momentum at values ranging from to 54. Inset: zoom- 
in of the spectrum with k/kp < 2. The dashed line indicate 
the ra value at which the effective mass at the Fermi surface 
shows divergence ~ 50. 



An obvious common feature of both 2D and 3D quasi- 
particle energy dispersion is that at some large mo- 
mentum, the dispersion shows a kink, which indicates 
the plasmon-resonance. This is the threshold wavevec- 
tor at which the quasiparticle can emit plasmons obey- 
ing energy-momentum conservation^^. This plasmon- 
emission- induced kink in i?k for k > kp occurs at all rs, 
with the threshold wavevector kth being rg-dependent. 
As Ts increases, this resonance feature becomes stronger. 
To further substantiate the connection between the dis- 
persion (and the associated effective mass divergence) 
and the plasmon resonance we show in Figs. [3 (2D) 
and IHl (3D) our calculated quasiparticle energy disper- 
sion together with the calculated quasiparticle damping 
Fk = /mE(k,^k)- It is obvious from Figs. [7| and |S1 that 
the "dip" (at k > kp) in the quasiparticle dispersion oc- 
curs precisely at the plasmon emission threshold wavevec- 
tor^^. The energy-momentum conservation makes it im- 
possiblei^ for plasmon emission by quasiparticles to hap- 
pen below the threshold momentum kth- Since the real 
and imaginary parts of the electron self-energy function 
are connected through the Kramers-Kronig causality re- 
lations, the plasmon emission threshold in the quasipar- 
ticle damping (i.e. ImS) shows up as a kink in the quasi- 
particle dispersion (i.e. ReE). This kink or the dip in the 
dispersion indicates a nonlinear instability of the quasi- 
particles over a finite range of momentum since the quasi- 
particle velocity becomes manifestly negative for k < kth 
where the dispersion is locally inverted. With increasing 
Ts, interaction effects become stronger and the plasmon- 
emission-induced kink also becomes deeper indicating a 
progressively stronger (with increasing rs) nonlinear in- 
stability around kth- The plasmon emission threshold 
kth{> kp) IS given by the solution of the following tran- 
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scendental equation: 



max 

kF<q<kt 



e(fc 



th 



2m 2m 



) = 0. 



(24) 



Below we discuss the (partial) connection between the 
plasmon emission phenomenon (leading to the 'kink' in 
the quasiparticle dispersion at the plasmon resonance 
momentum fc^^) and the dispersion instability, which is 
quite apparent in Figs. |S1 El [7| |H1 and at the same time 
emphasize some subtle features about this connection, 
which point to the qualitative nature of this connection. 
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FIG. 7: (Color online) 2D real and imaginary quasiparticle 
energy as a function of momentum at values ranging from 
2 to 20. 

A closer examination of the energy dispersion (see 
Figs. |S1 and |BJ provides interesting insights. Here we fo- 
cus on the momentum region where the dispersion shows 
instability. The divergence of the effective mass (i.e. the 
inverse slope of the energy spectrum) for quasiparticles 
at different wavevecter (not just at the Fermi surface) is a 
good indication of these instabilities (also see Figs. Eland 
[TUIl . Due to the presence of the plasmon resonance which 
we described in the last paragraph, one effective mass 
divergence happens at momentum (which we call the up- 
per instability momentum fc^) slightly smaller than the 
threshold resonance momentum (kth) for both 2D and 
3D systems, i.e. m*{k) ^ oo as k —> k^. The energy 
dispersion reaches a local maximum when k — k^. As we 
mentioned, this feature exists for all values. However, 
for low electron density, this is not the only instability 
feature that occurs. In fact in both 2D and 3D, for Ts 
above a certain value, a quasihole instability emerges at 
a smaller momentum, which we call the lower instability 
momentum fc^* . The energy dispersion reaches local min- 
imum when k = ki , and m*{k) —^ooask—^k*. In other 
words, the slope of the energy dispersion is negative in 
the region k < k^ , positive for k'^ < k < k^^, negative for 
k^ < k < kth, and then positive for k > kth- So it is 
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FIG. 8: (Color online) 3D real and imaginary quasiparticle 
energy as a function of momentum at Vs values ranging from 
3 to 54. 



clear that, for big enough (but still smaller than r*), 
there are two instability regions k < k^ and k'^ < k < kth 
away from Fermi surface (fc* < kp < k^), where the mo- 
mentum dependent effective mass is negative. (For very 
small Vs we only have the quasiparticle instability region 
k'^ < k < kth-) We note that while the quasiparticle 
dispersion instability at fc^ < kth is induced by plasmon 



emission due to the obvious connection between fc^ and 
kth, no such simple explanation seems to apply to the 
corresponding quasihole instability at A; < k^ . 

For Ts < r*, the lower and upper instability regions 
remain below and above Fermi surface, and therefore are 
not of any particular significance since the quasiparticle 
damping (Fk = ImE(k, ^k)) is finite. As increases, 
both of these two regions grow larger: fc^* increases, fc^ 
decreases, and kth increases with increasing r^. If one (or 
both) of these two instability regions reaches the Fermi 
surface, the effective mass on the Fermi surface diverges. 
This is indeed what happens except that the details are 
slightly different depending on the system dimensionality 
(2D or 3D). In 2D, as r^. — > r*, both instability regions 
reach fcj;^ at the same time, whereas in 3D it appears that 
the quasihole instability reaches kp first at — r*. The 
difference between 2D and 3D results arises from the very 
different density of states in the two cases. 



10, 

kf 



kf starts out at zero momentum. 



increases and at the same time 



For 2D, when ~ 
and as increases, 
fc^ decreases. As rg reaches the critical value r* ~ 18, 
ki = k'^ — kp^ which is the critical value at which the 
effective mass at the Fermi surface diverges. When > 
r*, there is no effective mass divergence in the whole 
momentum space, and the dispersion is inverted for all 
k < kth- The 3D case is a little different. As ~ 28, k'l 
first emerges at zero momentum, and keeps increasing as 
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rg increases while decreases. As reaches the value 



50, kr 



kf, and the effective mass at the Fermi 



surface diverges. However at this Vg, another effective 
mass divergence happens at k^ > kp- Only when 
increases to more than 58, k* and k^ meet each other 
at around l.Gkp, after which the effective mass does not 
diverge in the whole momentum space. 

Because it is difficult to visually locate the local max- 
imum and minimum in the energy dispersion, we show 
in Figs. El and II Ul our calculated momentum-dependent 
quasiparticle effective mass m*(k) — [k~^dE]f^/dk]^^ . Of 
course, the concept of a quasiparticle effective mass far 
away from the Fermi surface is not particularly mean- 
ingful since these quasiparticles are necessarily highly 
damped with short lifetimes. Nevertheless it is impor- 
tant to realize that the effective mass divergence initially 
develops for < r* (i.e. at densities much larger than 
the critical density for the effective mass divergence at 
the Fermi surface) at wavevectors k^> kp and /cj* < kp 
far away from the Fermi surface. Eventually at — r* 
the effective mass divergence reaches kp, becoming at 
the same time a dispersion instability where the quasi- 
particle dispersion essentially becomes almost flat all 
around kp, implying rn(k) to be divergent at the Fermi 
surface. This approximate "band-flattening" , i.e. ap- 
proximately independent of k, for large enough Ts values 
is also apparent in Figs. El El |H1 
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FIG. 9: (Color online) Momentum dependent 2D quasiparti- 
cle effective mass for Ts = 8, 12, 16, 20. 



FIG. 10: (Color online) Momentum dependent 3D quasipar- 



believe that, while the basic mass divergence and band 
flattening phenomena are generic phenomena in strongly 
interacting electron liquids, the precise value of r* must 
depend on the approximation scheme involved, and it is 
likely that the "correct" r* is larger than that obtained 
within RPA. Our reason for using spinless electrons in 
our calculation is that within RPA the electron system 
undergoes a ferromagnetic transitioriii of either Stoner or 
Bloch type to full spin polarization at a critical value 
which is lower than r* , the critical Vg for the mass diver- 
gence, and therefore it is more appropriate to carry out 
the dispersion-instability/mass-divergence calculation for 
a fully spin-polarized electron system as we do in the 
current work. We also notice that previous theoreti- 
cal worksiSiiS predicted possible partial spin-polarized 
ground state in 3D electron systems. We do not consider 
this effect in our calculation because within the same ap- 
proximation scheme as in the present work, we find^^ 
that in 3D electron systems, the densities of the partial 
spin-polarization region are much higher than the den- 
sity that corresponds to the divergence of effective mass, 
and therefore partial spin-polarization has no effect on 
the dispersion instability. Partial spin-polarization does 
not occur in 2D systems^®, and thus this issue does not 
arise for our 2D calculations. 



Before concluding this section we mention that the crit- 
ical r* for the effective mass divergence (18 for 2D and 
50 for 3D) obtained in this work is slightly larger than 
that obtained in Ref. 1 (16 in 2D and 48 in 3D) since 
the current work considers spinlessfor equivalently, spin- 
polarized) electrons whereas Ref. |l| dealt with a param- 
agnetic electron system with a spin degeneracy of 2. It is 
interesting that r* shows only a weak dependence on the 
spin-polarization properties of the electron liquid. We 



IV. DISCUSSION AND CONCLUSION 

We show in this paper that the standard (and widely 
used) ring diagram approximation (RPA) self-energy cal- 
culation leads to a quasiparticle dispersion instability in 
the strongly interacting electron liquid (both 2D and 3D) 
at a low critical density - the quasiparticle dispersion be- 
comes almost "flat" around kp for — r* and beyond 
the critical point, for > r*, the quasiparticle effec- 
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tive velocity becomes negative implying an instability. 
The band flattening phenomenon is similar to that envi- 
sioned by Khodel and ShaginyanS. in a different context 
some years ago. We have also identified a physical mech- 
anism leading to the dispersion instability and the asso- 
ciated effective mass divergence. We find that it arises 
partially from plasmon (or collective mode) emission by 
the quasiparticles - in some sense, the "recoil" associated 
with strong plasmon emission slows down the quasipar- 
ticles, eventually effectively "stopping" it, leading to the 
mass divergence and band flattening. We emphasize how- 
ever, that the plasmon emission and recoil mechanism is 
at best a partial physical mechanism for the dispersion 
instability - it plays a crucial role, but is not the whole 
story. 

Two important questions immediately arise in the con- 
text of our theoretical findings: (1) What does the dis- 
persion instability signify or imply? (2) What, if any, 
are its experimental implications and consequences? We 
cannot answer either of these questions definitively. But 
we can speculate some possibilities. Since the answer to 
the second question obviously depends on the answer to 
the first question, we first discuss some possible answers 
to the first question. 

The meaning or implication of the dispersion instabil- 
ity associated with the effective mass divergence and the 
quasiparticle band flattening has earlier been discussed 
in the literature. Two of us have earlier speculated^ that 
the effective mass divergence is essentially a continuum 
version of the Mott transition as envisioned, for example, 
in the Brinkman-Rice scenarioS&. Such a Brinkman-Rice 
Mott transition scenario had earlier been invoked^^ in the 
context of the "almost localization" phenomenon in nor- 
mal He-3 where the short-ranged inter-fermion interac- 
tion (in contrast to the long-ranged Coulomb interaction 
we are considering) is known to lead to quasiparticle mass 
divergence in the strongly interacting regime (which cu- 
riously happens at very high, rather than very low, den- 
sities in He-3 since the interaction is short-ranged in that 
system). Usually the continuum analog of the Mott tran- 
sition in an electron liquid system is thought to be the 
Wigner crystallization transition, which, according to the 
best current quantum Monte Carlo (QMC) simulations, 
occur at a critical of 38 (2D) and 90 (3D). It is entirely 
possible that the effective mass divergence and the band 
flattening we find within RPA is the precursor to (or per- 
haps the RPA signature for) the Wigner transition, but 
we cannot prove or establish this speculation with any 
kind of theoretical arguments. It would be highly desir- 
able in this context to carry out QMC calculations to 
study the dispersion instability, but QMC studies are re- 
liable only for the true ground state properties, and may 
be quite inaccurate for effective mass calculations. 

A second possibility, closely related to the Wigner crys- 
tallization transition discussed above, is that the dis- 
persion instability is essentially a charge density wave 
instability^*^. This is not an absurd proposition given 
that the possibility of a density-driven charge density 



wave instability in a jellium electron liquid was origi- 
nally discussed^ by Overhauser more than forty years 
ago. We have therefore carefully analyzed the wavevector 
dependent contributions to the quasiparticle self-energy 
as well as calculated the Fermi liquid interaction function 
/(k, k') to see if there are any characteristic wavevectors 
which predominantly contribute to the dispersion insta- 
bility. As should be obvious from our results, there are 
no characteristic wavevectors in the effective mass diver- 
gence phenomenon that we have discovered since it arises 
from a dispersion instability in which the whole disper- 
sion is strongly affected and all wavevectors seem to be 
equivalent. The plasmon emission threshold wavevector 
kth seems to be special since the "seed" or "source" of 
our dispersion instability at least partially lies in the soft 
plasmon emission process, but we see no reason to as- 
sociate kth with a characteristic charge density wave in- 
stability. Similarly, the upper (fc^) and the lower (fc^*) 
critical wavevectors where the mass divergence first man- 
ifest itself (Figs. El and EI) could be characteristic charge 
density wave vectors but we have seen no theoretical in- 
dication of such a CDW instability, at least within our 
RPA theory 

Anther possibility that has been much discussed in the 
literature 2.4,5,6 ^ superconducting or fermionic con- 
densation instability associated with the quasiparticle en- 
ergy degeneracy in the flat dispersion around kp. Such 
a superconducting instability is entirely different from 
the usual Kohn-Luttinger superconducting instability (in 
some high angular momentum channel) that is known to 
exist (with an exponentially low superconducting transi- 
tion temperature) in interacting electron liquids. In our 
problem of dynamically screened Coulomb interaction, 
such a superconducting instability could possibly arise 
from the exchange of virtual plasmons (i.e. plasmon- 
mediated superconductivity) since emission of real plas- 
mons is an underlying mechanism for our dispersion in- 
stability. We have, in fact, little to add to the existing 
discussion in the literatureSi^i^ on the issue of a super- 
conducting instability in the context of the dispersion 
instability except to note that the transition tempera- 
ture is likely to be rather low for such a superconducting 
state. 

We note in this context that Nozieres carried outi^ a 
penetrating and trenchant analysis of the original band- 
flattening phenomenon introduced in Ref. U While dis- 
cussing the instability, Nozieres^ was quite pessimistic 
about the theoretical possibility of such an instability 
existing in a realistic model that allows for screening. In 
fact, Nozieres concluded that "screening of a strong long 
range interaction is such that the instability threshold 
cannot be reached". Our work shows that this conclu- 
sion of Ref. 5 is, in fact, too pessimistic since our one-loop 
calculation is precisely an expansion in the dynamically 
screened long-range Coulomb interaction. Thus, the ef- 
fective mass divergence can certainly occur even when 
screening of a strong long range interaction is explicitly 
incorporated in the theory in contrast to Nozieres' con- 
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elusion in Ref . H. 

A real theoretical concern is the possibility that the 
effective mass divergence (and the dispersion instability) 
is just an unfortunate artifact of our specific self-energy 
approximation, namely the single-loop (leading-order in 
dynamically screened interaction) RPA self-energy calcu- 
lation. Although such a possibility can never be ruled out 
theoretically (short of an exact treatment of the strongly 
interacting electron liquid problem) we have arguedi else- 
where that this is unlikely to be the case here, i.e., there 
is very good reason to believe that the dispersion in- 
stability at a low electron density is a generic property 
of the strongly interacting electron liquid (but the ac- 
tual value of Ts depends on the approximation scheme). 
Without repeating the arguments already discussed in 
Ref. Jj, we point out that RPA (which is a self-consistent 
field approximation, not a perturbative expansion in 
although RPA does become exact in the high density 
limit) works well (as compared with experimen- 
tal data) at metallic densities {rg ~ 3 — 6) in 3D systems^ 
and at very low densities (r^ '--^ 10 — 20) in 2D electron 
systemsiS^,. Second, the RPA prediction for a ferromag- 
netic instability in electron liquids turns out to be gener- 
ically "correct" , i.e. the most accurate QMC calculations 
also predict ferromagnetic instabilities in 2D and 3D elec- 
tron systems except that RPA underestimates the critical 
Ts for the electron liquid ferromagnetic instability. This 
suggests that the RPA prediction for the dispersion in- 
stability is also likely to be generically correct. Third, 
the original dispersion instability, envisioned in Ref. 2, 
used a toy model (involving a rather unrealistic long- 
range interaction) which is completely and fundamen- 
tally different from our dynamically screened Coulomb 
interaction model - the fact that two completely different 
interaction models arrive at very similar qualitative con- 
clusions about the dispersion instability (both in 2D and 
3D) is again strongly suggestive of the possible generic 
nature of the instability. Fourth (and perhaps the most 
important of all), there are no general theoretical argu- 
ments against such a dispersion instability, and there- 
fore, if the interaction is strong enough (i.e. large 
enough), there is no theoretical reason for the instability 
not to occur. This point becomes even stronger in the 
context of the existence of such an effective mass diver- 
gence ("almost-localized Fermi liquid") in normal He-3 
interacting via the short-range interaction. We also note 
that, as shown in Ref.Q, the quasiparticle effective mass 
divergence occurs in many-body approximations going 
beyond the RPA (e.g. the Hubbard approximation^'^), 
and the same is obviously true for the dispersion insta- 
bilities. Thus the instability exists in self-energy calcula- 



tions beyond the one-loop approximation. 

Finally, we discuss the experimental implications of our 
results. First we note that the critical density (r* — 18 
(2D); 50 (3D)) involved in the dispersion instability is 
extremely low for meaningful experiments to be carried 
out in real systems. Also, our T = theory does not 
say much about the finite temperature situation where 
these experiments are necessarily carried out. In fact, 
the temperature scale at such low densities is likely to be 
small, and therefore one may have to go to unrealistically 
low temperatures to see any experimental consequences 
of the dispersion instability even if the instability is real. 
It is important to emphasize that our RPA values for 
critical r* are probably lower bounds on the true r* which 
is hkely to be larger than r* = 18 (2D) and 50 (3D). 

With these caveats in mind, it is interesting to note 
that there have been several recent experimental claims 
of the observations^ of quasiparticle effective mass diver- 
gence in Si MOSFET-based low-density interacting 2D 
electron system at ^ 10. These claims are, however, 
quite controversial, and in more dilute 2D hole systems, 
where could reach as low as 20—30, no such mass diver- 
gence has been reported. We believe that these recently 
reported effective mass divergence claims in Si MOSFETs 
are most likely not connected with the dispersion insta- 
bility discussed in our work. This is particularly true in 
view of the fact that the reported effective mass diver- 
gence does not seem to be a dispersion instability, and 
furthermore, the experimental critical r* (~ 10) is far 
too low compared with the theoretical finding 18). In 
fact, we expect our theoretical r* to be much lower than 
the RPA value (r* ~ 18 in 2D) predicted in our work once 
quasi-2D finite width effects and many-body effects go- 
ing beyond RPA are included in the theory. In addition, 
the semiclassical experimental technique of using Dingle 
fits to temperature-induced SdH amplitude decay^** in 
extracting the quasiparticle effective mass is highly sus- 
pect in such a strongly interacting quantum system. The 
experiments claiming the effective mass divergence also 
completely ignored the strong temperature dependences^ 
of the quasiparticle effective mass which certainly invali- 
dates the semiclassical Dingle fitting procedure employed 
in the experiments in obtaining the quasiparticle effec- 
tive mass. We therefore conclude that the experimental 
consequences of the dispersion instability in low-density 
interacting electron liquids remain an interesting open 
challenge for the future. 
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